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Abstract 

The nuclear reaction e+p + d — > ^He + e is considered at thermonuclear ener- 
gies. The motion of the electron is treated within the adiabatic approximation 
and the pd scattering state is constructed in the form of an antisymmetrized 
product of the bound state wave function of the deuteron and of the wave 
function of the pd relative motion. The latter is calculated using an effective 
pd potential constructed via the Marchenko inverse scattering method. The 
bound state wave function of ^He is obtained using Faddeev-type integrod- 
ifferential equations. The reaction rate thus obtained for the solar interior 
conditions is approximately 10""^ of the corresponding rate for the radiative 
capture pd — > ^He7. 
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I. INTRODUCTION 



Burning of hydrogen in the main sequence stars mainly occurs through the pp-chain 
which begins with the reaction pp — > e^dv. It is generally accepted that the second 
step of this chain is the radiative capture of protons by deuterons, 

p + d — >^He + 7. (1) 

However, due to the high densities in stars, the helium nuclei and other intermediate products 
of the pp-chain can emerge not only from two-body but from three-body initial states as 
well. Thus, for example, ^He can be formed in the radiative capture (|l|) as well as in the 
three-body nonradiative fusion process 

e + p + d — > ^He + e . (2) 

The knowledge of the production rates of various nuclei in the stellar plasma is very im- 
portant, not only in understanding the production of energy in stars but also in explaining 
the abundance of the elements observed in nature, and in describing events during the first 
thousand seconds of the evolution of the universe which were predominantly determined by 
the nucleosynthesis of light elements [@]. The abundance of light elements in the Universe, 
together with the Hubble expansion and the relic backround radiation are experimental ev- 
idences for the idea of the hot origin of the Universe in the Big Bang process which in turn 
is closely related to the Grand Unification and QCD theories 0,^. 

Values of the primordial abundance of light elements put some constraints on the baryon 
density of the Universe as well as on the number of species of light particles. For instance, 
present theoretical estimates for the primordial abundance of d, ^He, ^He, and ^Li can be in 
agreement with the corresponding experimental values only if the number of the neutrino 
species is < 3.9 0]. Furthermore, from the theory of nucleosynthesis one can derive a 
stringent limit to the existence of new light particles and even a bound to the mass of the 
r-neutrino, namely, between 0.5 MeV and 30 MeV 0J§]. 

Any theory on the evolution of the universe or nucleosynthesis must deal with the total 
rates of nuclei production. In this respect the thermonuclear reactions with two-body initial 
states, such as the reaction (||) have been extensively investigated 0. However, the role of 
the three-body mechanism in nucleosynthesis has not yet been properly studied despite the 
fact that the three-body processes have different selection rules and due to this it can have 
an influence on the production of light nuclei. 

The aim of this work is to estimate the relative significance of the nonradiative process 
(0) in comparison with the radiative capture (|I]) in the stellar pp-chain. In stellar plasma 
nuclei are surrounded by electron gas which has a twofold influence on nuclear fusion pro- 
cesses. We can distinguish between static and dynamic electronic effects [Q. The former is 
the screening of the Coulomb repulsion between the nuclei, and the latter stems from the 
coupling between the electronic and nucleonic degrees of freedom. Due to this coupling, en- 
ergy and angular momentum can be transferred from the nucleons to the electrons according 
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to the prevailing conservation laws. These dynamic electronic effects can pave the way for 
a variety of possible fusion reactions in the plasma, which otherwise are forbidden. 



The nonradiative fusion (H) is an example where these effects are manifested. The elec- 
tron which is in the vicinity of the pci-pair while they are interacting can carry away the 
excess energy, leaving the three nucleons in a bound state. This is a kind of Auger transition 
in the continuous spectrum. Since electrons move much faster than nucleons, in considering 
the scattering of an electron by a pd-pair, its motion can be treated within the adiabatic 
approximation. In this approximation the nucleons are considered as being fixed at their 
spatial positions during the electron scattering and thus the corresponding amplitude will 
depend on the nucleon coordinates parametrically. The physical amplitude can then be 
obtained by averaging over these coordinates with the help of the wave functions describing 
the initial and final configurations of the nucleons. To obtain the amplitude of the reaction 
(0) we shall average the fixed-scatterer amplitude over the pd scattering state and the ^He 
bound state. 

The paper is organized as follows. In Sec. II we describe our formalism and outline the 
procedure employed to evaluate the various ingredients used to obtain the reaction rate. In 
Sec. Ill we present our results and conclusions. Some details concerning the derivation of 
the reaction rate formula are given in the Appendix. 



We are concerned with electron collisions with nuclei in a stellar plasma consisting of 
protons, deuterons, and electrons. Let \l/k be the wave function of the relative motion of the 
pd-pa.iT with momentum k and let p be the momentum of the electron with respect to this 



represent particle densities. 

In stellar plasmas the momenta k and p are distributed according to Maxwell's law 



where A^k and A^p are the probability densities, is the proton-deuteron reduced mass, m 
is the electron mass, k is the Boltzmann constant, and G is the plasma temperature. Thus 



II. FORMALISM 




where the states in the continuum are normalized as 
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the reaction rate must be averaged over the initial momenta k and p and integrated over 
the final momentum p', i.e., 

{'^) = J J Jdkdpdp'n{k,p^p')N^Np. (4) 

In what follows we shall discuss the various ingredients of this formula in somewhat more 
detail. 



A. Transition operator 

Consider a four-body system consisting of an electron and three nucleons described by 
the Jacobi vectors shown in Fig. I. The total Hamiltonian of this system consists of three 
terms, namely, 

H = Hs + ho + Ve, (5) 

where H3 is the total nuclear Hamiltonian, ho is the kinetic energy operator for the free 
motion of the electron with respect to the center of mass of the nucleons, and Ve is the sum 
of the Coulomb potentials between the electron and nucleons 

T/ ^1^' ^26^ Z,e^ 

Ve = . (6) 

n r2 rg 

Here r, are the distances between the electron and the nucleons, and Zi is the charge-operator 
for the i-th nucleon with 

where |l/2, Tj > is the isotopic state of the i-th nucleon with Xj being the third component 
of its isospin. 

Since the motion of nucleons is much slower than that of the electron, the T-matrix 
describing the electron scattering can be found with the help of the fixed-scatterer approx- 
imation 

T{z) = Ve + VeGo{z)T{z) , (8) 
where z is the total energy and 

Goiz) = {z- ho)-' , 

is the free Greens' function. Parametrical dependence of the potential @ on the nuclear 
Jacobi vectors {p, r} and on the nuclear isospin state \ri > of the three nucleons, makes the 
fixed-scatterer T-matrix (|^) also parametrically dependent on them, 

(p', r'; T]'; p'\T{z)\p, r; rj; p) = 6.^^'6{p' - p)6{r' - r)T^,p(p, r; z) , 
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where p and p' are the initial and final momenta of the electron. The physical T-matrix is 
obtained by the averaging, 

(^3,plT(z)|v[/k,p) = E / / ^P^r^3"'^*(p,r)Tp%(p,r;z)^f (p,r) , (9) 

where ipl^ and ip^^ are the spatial parts of the three-nucleon wave functions corresponding 
to the spin-isospin state Ix'? >, 

\XV >= |((siS2)si2, S3)^(r; ((tit2)ti2, ^s)^^) (10) 



(p,r) = (p,r;xr?|^k), (11) 

and 

^3^^(p,r) = (p,r;xr7|*3). (12) 



B. Proton deuteron scattering state 

The thermonuclear energies we are considering are far below the threshold of the deuteron 
break-up and we can therefore neglect the three-nucleon continuous spectrum (1 + 1 + 1) 
and construct the (1 + 2) scattering state in the form of the antisymmetrized product of the 
proton-deuteron relative motion wave function and the deuteron wave function iIj2, 

*k = A{(pkip2XV} ■ 

Since the deuteron wave function is antisymmetric, the antisymmetrizer A involves only the 
permutations P23 and P13, 



A 



V3 



'1 - P23 - 



13) 



(13) 



Both (y9k and ip2 are assumed to have S-wave components only with the total spin of the 
nd-system equal to 1/2. Before the antisymmetrization, the spin-isospin states of the nd- 
system are \xv) = |(('5iS2)ls3)|cr; ((tit2)0t3)|r). Thus after the nucleon permutations we 
obtain 



1 

V3 



¥'k(p)^2(r) 



(14) 



- <^k(--p + ^r)^/'2(p + -r) 

- V?k(--P- ^r)^/'2(-r-p) 



((siS3)ls2)^a;((M3)0t2)^r 
((s3S2)lsi)^(T;((t3t2)0ti)^r 



The spatial components i/j^^{p,r), Eq. (pTj), can be obtained by projecting on to the spin- 
isospin states ([T0|). 
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To find the deuteron wave function '?/'2(i") we solve tlie two-body Sclirodinger equation 
with the Malfliet-Tjon I-III A^A^-potential. In order to obtain the relative motion wave 
function v^klp) the two-body scattering problem is solved using the Jost function method 



proposed in Ref. |Tl[] with an effective proton-deuteron potential Vpd{p) which consists of 
two terms: 



vM = v;(p)exp -— + Vsip). (15) 




The first term describes the proton-deuteron Coulomb interaction. Since deuteron is not a 
point-like particle, we take into account the spherically symmetric distribution of its charge 
in space by using 



Vcip) = / dr [rilj2{r)] 

p Jo 

At large p (beyond the deuteron radius), Vc{p) coincides with the simple Coulomb potential 
e^/p. However at small p it behaves quite differently. In particular, Vc{p) is not singular 
at p = 0, but instead V^(0) = 0. The exponential factor stems from electron screening and 
follows from the standard Debye-Hiickel theory |jl2|]. For the Debye radius we use the value 



Pd = 21800 fm 

which corresponds to the solar plasma conditions and is typical for other stars [§]. 

The Vs{p) is the strong (nuclear) pd interaction and we constructed it using the i- de- 
pendent Marchenko inverse scattering method ll^Jl^ which we briefly describe next. 



C. Marchenko Inverse Scattering method 

In the Marchenko inverse scattering method a unique, energy-independent, ^-dependent, 
local potential Vf {p) can be constructed which is phase equivalent to the nd doublet channel 
effective interaction at all energies. This potential is obtained from 

V,ip) = -2^^K,ip, p) (16) 

where the kernel Ki{p,p') obeys the Marchenko fundamental equation 

Ke{p, p') + F,(p, p') + / K,{p, p")F,{p\ p')dp" = 0. (17) 

Here, p is the relative distance between the neutron and deuteron and is canonically conju- 
gate to the relative momentum k. The driving term Fi{p^p') is given by 

1 /• + 00 

Fi{p. P) = 7r <(^P)[l - S^{k)]wt{kp')dk + A^wiMwiM. (18) 

Si{k) is the S-matrix for the specific partial wave and the function wl^z) is related to 
the spherical Hankel function h^i^\z) by 
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wi{z) =ie'''^zhi'^\z). (19) 

Furthermore, Ai is the so-called asymptotic bound state normalisation constant, while 
-4M£'^/3, being the three-body bound state energy and M the nucleon mass. 

As is apparent from its definition (p!8| ) the evaluation of Fe{p,p') requires the knowl- 
edge of the S-matrix for all real energies from the elastic scattering threshold to infinity, 
together with the refiection property Si{—k) = 1/Si{k), as well as the binding energy and 
the corresponding asymptotic bound state normalisation constant. It is greatly simplified 
by choosing a rational (Bargmann-type) parametrisation 

s,(k) = ^ n ^ (20) 

The number A^^ is to be taken odd to satisfy the requirement Se —>■ 1 for k 0. The are 
complex numbers used to fit the (numerically) given S-matrix. 

With the choice (^) the integration in Eq. ( [TE| ) can be easily performed analytically, 
the result being 

F,(p, p') = -iY. «(«mPX(«^p') + A^wi{htp)wt{htp') . (21) 

m=l 

Here, i?^ are the coefficients of the residues of the integrand for the poles of ( pOD which 
are lying in the upper half complex k-plane (excluding the one corresponding to the bound 
state at = ib^). The separable form of Fe{p, p') provides us with an algebraic solution of 
the integral equation (|l^ for the kernel Ki{p,p') from which the potential can be obtained 
viaEq. (|T|). 

In order that the potential be unique, the asymptotic bound state normalisation constant 
Ai must be given a definite value. We choose it to be 



fpi—ihp] , . 



where fi is the Jost function, and fi{k) = dfi{k)/dk. 



D. Wave function of ^He 

The integrodifferential equation approach to few- and many body systems developed by 
Fabre de la Ripelle and collaborators is used to construct the bound state wave 

function of ^He. In this method the potential U is written as a sum of two-body interactions 

f^ = Enr.,), (23) 

i<j 

and the three-body bound state wave function is written as a sum of two-body amplitudes 
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i<j 



(24) 



where i?[Lm] (^) is the harmonic polynomial of lowest degree [L^] occuring in the harmonic 
polynomial expansion of the wave function, and x represents the nucleon coordinates Xj with 
Tij = |xj — Xjl, and tq = [2/3Y,i<j ^ij]^''^ hyperradius. 

With the above expansions one has to solve, instead of the Schrodinger equation, the 
Faddeev-type equation for the amplitude F{rij,ro) 



{T-E)Hu.{^)F{r,^,ro 



-V{r,,)H[L^]{^)Y.F{ne,ro) 



(25) 



A solution of this equation will be an approximate solution of the Schrodinger equation for 
two-body amplitudes where pairs are in ^-states. Another Faddeev-type equation can be 
obtained by extracting the hypercentral potential ViL,„](ro) of V{rij) [0] and by writing 

A{A- 



(T + 



-V[L,„](ro)-i?)%„](x)F(r,„ro) 



(26) 



-[^(%) - V[L^]{ro)]H[L^]{^)Y.F{ne,ro) . 

k<£ 



This equation forms the basis of the integrodifferential equation approach (IDEA) to few- 
and many-body systems |jl6| and takes into account in an approximate way, via the hy- 
percentral potential V[L,r,]{^o), the effects of the coupling between the orbitals £ 7^ of the 
spectator particle and the interacting pair. For [L^] = then if[£,^](x) = 1. We notice that 
by summing over all pairs one generates the Schrodinger equation but here the two-body 
potential is the residual interaction on the right hand side of (p6|). 



We assume that we have a central spin-dependent nucleon-nucleon potential of the form 



where the projection operators P^^ and Pf^ are acting on the singlet- and triplet-even 
states respectively. In order to proceed two further steps are required. In the first step 
F{rij,ro) is written, 

F(r,„ro) = P(0„ro)/rf-^)/^ 

where Qj = 2rfj/rQ — 1 and D = 3A — 3. In the second step Eq. (^) is projected on the r, 
space to give two coupled integrodifferential equations 



— ^ 



m 



-Vo{ro) + E 



PoiC^ro) 



+ 



rfe2 



m 



^0 



AiA-l) 



Vo{ro) + E 



^o"'(C,ro) 



+ 









2 

- -1/3+ 




2 




+ V^+ 




2 
- 



- Voiro) 
<(C,ro 



no^(C,ro) 



(27) 



- Voiro) 

n^(Cro) 



no"'(C,ro) 
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where 



n^(C,ro) = Po"(C,ro) + J ^ /("o)(C,C')^o"(C',ro)rfC' (28) 

with n = S, S'. 

The Vq is given by 

^^^^TT^ + ^oWoi^^ ~ ^'4' ^^^^ 

with Cq = {D — 3)/2, while the weight function Wq by 

w^o(C) = (i-C)^^-')/'(i + C)'/'. (30) 

The \/« = V» (roy^(l + C)/2), t = 1+, 3+, are the even singlet and triplet nucleon- 
nucleon potentials. The kernels f(})){(,(')y n = S,S' result from the projection on to the 
Tij space. More details concerning these kernels and other technical points can be found in 
Refs. IT6IJTT 



Once the components Pij{C,fo) are found, one can construct the various symmetries for 
the bound state wave function \&3 ( [^]). However, in the present calculation we used only 
the space symmetric part i.e.. 



(Ariz's) = 



V2 



It is emphasised here that the IDEA for the three-body case and for S-projected potentials, 
is equivalent to the exact Faddeev equations. 



E. Electron scattering 

The scattering of the electron on the i-th nucleon, is described by the t-matrix 

U{z)=V,^V,G^{z%{z), (31) 

where 

and Gq{z) is the free Green function. Using ( pi]) we can rewrite Eq. (||) as a sum of Faddeev 
components, viz., 

3 

T = Y.% (32) 

1=1 



with 



T, = t, + t,GoE^i (33) 



9 



Iterative solution of these equations provides the multiple scattering series 

T = ti+t2 + h + hGot2 + hGoh + tiGots + hGoh + hG^h + hGot2 + ■■■ . (34) 

Since the average energy of the Coulomb interaction of the electron with the nucleons is 
of atomic order of magnitude (~ 10 eV) and the average collision energy is {kQ) ~ 10^ eV, 
we can omit the higher order rescattering terms in Eq. (PW), that is, we use the following 
approximation |1^ 



T ^ti+t2 + h. (35) 
For the same reason we can apply the Born approximation for the Coulomb t-matrix, i.e. 

Qe2 



ic(p',p;<5) 



27r2 (p - p'; 



,2 ' 



where tc is the two-body Coulomb t-matrix for an electron scattered off the charge Qe. 
To obtain matrix elements of the T-operator (B^), we use the following basis states 



|p,r,p; xv) 



p,r,p; ((siS2)si2S3)^cr; ((^1^2)^12^3)^^ / , 



(36) 



where p is the electron momentum with respect to the center of mass of the nucleons. Using 
the space displacement operators exp(— zAjp), where Aj is a vector directed from the 
center of mass to the i-th nucleon, we also construct the shifted basis states, 

|p,r,Pi;x^) =exp(-iAip) |p,r,p;xr/) , 

which differ from (|3^) in the sense that the electron has the same momentum Pi = p but 
with respect to the nucleon (not to the center of mass). Thus, we obtain 



with 



(p', r, p'; xV \T{z)\p, r, p; xv) = ^x'x^v'v^ip' - p)S{r' - r)Tl, Jp, r; z 



T^ip, r;z)=Y^ exp [z(p - p') A J t,(p', p; Ql) , 



(37) 



(38) 



m=l 



where the effective charges are 



Using the 6j-symbols |^ for recoupling the isospins, we obtain 



1 

2' 


Ql = 


1, 


5 


Ql = 


1 


6' 


3 



where the upper index designates the |((tit2)0t3)^^) and 1 the \{{tit2)lt^)ly states 
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III. RESULTS AND CONCLUSIONS 



The effective pd potential obtained using the Marchenko inverse scattering theory, de- 
scribed in Sec. Q for the £ = partial wave, is plotted in Fig. II. This potential reproduces 
the experimental nd phase shifts as the Coulomb pd interaction has been treated separately 
as described in Eq. ([T5|). Since these phase shifts, however, are available only at low energies 



we use for large values of k the phase shifts obtained via the Faddeev equations |2T[]. The 
oscillations in the interaction region are due to the opening of the break-up channel as well 
as to the behaviour of the phase shifts at large k values {k > 1000 MeV). 



The wave function of the ^He-nucleus was obtained by solving the system (|2^) with the 
Malfliet-Tjon I-III (MT I-III) nucleon-nucleon potential as input. The binding energy 
obtained is 8.86 MeV while the root mean square radius is 1.685 fm. However, for our final 
calculations we used the energy release in the reaction (|^), viz., 

£ = E3 — E2 , 

correponding to the difference between the experimental binding energies of the deuteron 
and ^He, namely E2 =2.224574 MeV and =7.718109 MeV fS^]. 



The rate of the reaction can be presented in the following factorized form (see the 
Appendix) 

(7^) = npUdTieiT,) , 

where the quantity (S) is analogous to {av) generally used in two-body reaction theories 
0]. The calculated rates for different temperatures of the plasma are presented in Table I, 
where the results are given in units of cm^mole~^sec~^ [^. These units are obtained when 
instead of UpUdne we multiply (S) by A^^, i.e. the Avogadro number squared. The second 
column of the Table I contains the values of 

7^e = iVl(E) 

while the third contains the rate for the radiative capture (|l]) which is also presented in the 
form 

= NA{av) , 

instead of npnd{av), that is in cm'^mole~^sec~^ units (the data are taken from Ref. [p3[] ). 

In order to assess the importance of the nonradiative process (|^) as compared to the 
radiative capture (|l]) when '^He nuclei are generated, we must compare the rates for these 
two processes per cm^ per sec. This requires the knowlege of the particle densities Up, na, 
and He, and therefore we need to specify the plasma conditions. However, by considering 
the ratio of the two rates, viz., 

ratio = He-. — r . 

{av) 

only one unknown parameter remains, namely the electron density He- We calculated this 
ratio for the reaction rates for (^) and (|I]) by using the value ng = IOOA^a cm~^ corresponding 
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to the solar interior plasma [|T| . The calculated ratio is given in the fourth column of Table 
I. The ratio for any other electron density Ue can be obtained by simply multiplying these 
values by ne/lOOA^A- 

It is seen that in the solar pp-chain the nonradiative fusion (Q) plays, apparently, a minor 
role. However, at the early stages of the universe when n^jNA ^ 100, this reaction must 
have been significant. 
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APPENDIX A: REACTION RATE: EXPLICIT FORMULA 



Due to the very low nuclear energies involved we shall use only S-wave components of 
the wave functions describing the deuteron, ip2, the pd relative motion, ip^, and the nucleus 
^He, ips- These functions depend only on the moduli of the corresponding vector variables. 
However, due to the antisymmetrization, the dependence on the angle between the vectors 
p and r is also present. 



The fixed scatterer T-matrix (|^) depends on all four angles that define the directions of 
the vectors p and r which are involved in (|37|) in the form of the following combinations 

Ai = -r + -p, 

^^ = 3^-2^' 
2 

A3 = -3P. 

Thus, the transition matrix (^) is a a 6-dimensional integral. The final averaging procedure 
involves 9 more integrations, 

{TZ) = npndn,{27,y J J J dkdpdp' 6{Ef - E,) \{^^,p'\T\'^^,p)\^ . 

The transition matrix {'^3,p'\T\'^y_,p) depends on (p' — p) and |k|. Moreover the Maxwell 
distributions are isotropic. Hence, the 5-dimensional integration over dQ^dQpH 



^ p is trivial 

and results in the factor 327r^. One more integration, over the final momentum, is performed 
with the help of the energy conserving 5-function and the 9-dimensional integral is reduced 
to a 3-dimensional one. Thus, finally we obtain 



(7^) = npUdUeiT.) 



where (S) is defined by 
(E) 



3\/mp^K^Q^ Jo 



dx ^ 



dkdp k'^p^ exp 

f-27r j-2-K 

/ di 
Jo 



p" 



drdp 

iPk (v'SrVie + pV4 - 3r.p/4 ) 



fi 



1 



dydzr p ^psip^r) 



+ ifk ( V9r-Vl6 + p2/4 + 3r-p/4 j ^2 

fiqpy 



■ (iqpy iqrz\ 



+ 3 exp 



r2/4 + p2 -|- Y-p 

r'^ /A + p2 — Y-p 

iqrz 



V 3 



+ <^fc(p)V'2(p) 



■ (iqpy iqrz 



+ 2 exp 



f npy 

V 3 



+ 2 exp 

iqrz 



-i2qpy\ 
3 J 

+ 4 exp 



-i2qpy\ 
3 J 



where the scalar product r -p is 



V -p = rp 



'^{l-y^){l-z^) 



Z^ cos 0r — 



+ yz 
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while 



and 



m . 



J2m£ + -k^+p\ 



q= jQ2+p2_2Qpx. 
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FIGURES 



FIG. 1. The Jacobi coordinates for the e + p + d system. 
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FIG. 2. The nd effective interaction obtained using the Marchenko inverse scattering method. 
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TABLES 



TABLE L Nonradiative, TZe = {Tl{pde — > e^He)), and radiative, TZj = {TZ{pd — > ^He7)), 
capture rates as functions of the plasma temperature Tq. The fourth column shows their ra- 
tio for Ue = lOOiVy^cm"^. The temperature Tg is in 10^ °K units, and the reaction rates in 
cm^mole~^sec~^ and cm^mole~^sec~^ respectively. The data for {lZ{pd — > ^He7)) are taken from 
Ref. |3|. 
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